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Abstract 

The entire function (of z) 

oo z k 

named after the British mathematician E.M. Wright, has appeared for the 
first time in the case p > in connection with his investigations in the 
asymptotic theory of partitions. Later on, it has found many other applica- 
tions, first of all, in the Mikusihski operational calculus and in the theory of 
integral transforms of Hankel type. Recently this function has appeared in 
papers related to partial differential equations of fractional order. Consider- 
ing the boundary-value problems for the fractional diffusion-wave equation, 
i.e., the linear partial integro-differential equation obtained from the classi- 
cal diffusion or wave equation by replacing the first- or second-order time 
derivative by a fractional derivative of order a with < a < 2, it was 



found that the corresponding Green functions can be represented in terms 
of the Wright function. Furthermore, extending the methods of Lie groups 
in partial differential equations to the partial differential equations of frac- 
tional order it was shown that some of the group-invariant solutions of these 
equations can be given in terms of the Wright and the generalized Wright 
functions. In this survey paper we consider some of the above mentioned 
applications of the Wright function with special emphasis of its key role in 
the partial differential equations of fractional order. 

We also give some analytical tools for working with this function. Be- 
ginning with the classical results of Wright about the asymptotics of this 
function, we present other properties, including its representations in terms 
of the special functions of the hypergeometric type and the Laplace trans- 
form pairs related to the Wright function. Finally, we discuss recent results 
about distribution of zeros of the Wright function, its order, type and indi- 
cator function, showing that this function is an entire function of completely 
regular growth for every p > — 1. 

Mathematics Subject Classification: primary 33E20, secondary 33C20, 30C15, 
30D15, 26A33, 45J05, 45K05 

Key Words and Phrases: Wright function, indicator function, asymptotics of 
zeros, generalized hypergeometric function, diffusion-wave equation, Green 
function, scale-invariant solutions, Erdelyi-Kober operators 

1. Introduction 

The purpose of this survey paper is to outline the fundamental role of the 
Wright function in partial differential equations of fractional order, to con- 
sider some other applications of this function and to give its analytical prop- 
erties including asymptotics and distribution of its zeros. Partial differential 
equations of fractional order (FPDE) are obtained by replacing some (or all) 
derivatives in partial differential equations by derivatives of fractional order 
(in Caputo, Riemann-Liouville or inverse Riesz potential sense). Mathe- 
matical aspects of the boundary-value problems for some FPDE have been 
treated in papers by several authors including Engler [5], Fujita [7], Goren- 
flo and Mainardi p], Mainardi [20]-[23], Podlubny [32], Priiss [31], Saichev 
and Zaslavsky [35], Samko et al. [36], Schneider and Wyss [37] and by Wyss 

From the other side, some FPDE were successfully used for modelling 
relevant physical processes (see, for example Giona and Roman [9], Hilfer 
[15] . Mainardi [21], Metzler et al. [26], Nigmatullin [28], Pipkin [31], Pod- 
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lubny [32] and references there). In applications, special types of solutions, 
which are invariant under some subgroup of the full symmetry group of the 
given equation (or for a system of equations) are especially important. 

Recently, the scale-invariant solutions for time-fractional diffusion-wave 
equation (with the fractional derivative in the Riemann-Liouville sense) and 
for the more general time- and space-fractional partial differential equation 
(with the Riemann-Liouville space-fractional derivative of order (3 < 2 in- 
stead of the second order space derivative) have been presented by Buckwar 
and Luchko pQ and by Luchko and Gorenflo [19], respectively. The case 
of the time-fractional diffusion-wave equation with the Caputo fractional 
derivative has been considered by Gorenflo, Luchko and Mainardi |13j . 

The plan of the paper is as follows. In Section 2, following the pa- 
pers by Djrbashian and Bagian [4], Gajic and Stankovic 8\, Luchko and 
Gorenflo [T9], Mainardi [23], Mainardi and Tomirotti 24J, Mikusihski [27], 
Pathak [30], Pollard [33], Stankovic [39], and Wright [42], [44], we recall 
the main properties of the Wright function including its integral representa- 
tions, asymptotics, representations in terms of the special functions of the 
hypergeometric type and the Laplace transform pairs related to the Wright 
function. Finally, we discuss new results about distribution of zeros of the 
Wright function, its order, type and indicator function, showing that this 
function is an entire function of completely regular growth for every p > — 1. 

In Section 3, we outline some applications of the Wright function, begin- 
ning with the results by Wright [H] in the asymptotic theory of partitions. 
Special attention is given to the key role of the Wright function in the theory 
of FPDE. Following Gorenflo, Mainardi and Srivastava |14j and Mainardi 
[20]-|23|. we consider in details the boundary- value problems of Cauchy and 
signalling type for the fractional diffusion-wave equation, showing that the 
corresponding Green functions can be represented in terms of the Wright 
function. We present also some results from Buckwar and Luchko [I], Goren- 
flo, Luchko and Mainardi [13J, Luchko and Gorenflo [19j concerning the ex- 
tension of the methods of the Lie groups in partial differential equations to 
FPDE. It will be shown that some of the group-invariant solutions of FPDE 
can be given in terms of the Wright and the generalized Wright functions. 

We remark finally that the present review is essentially based on our 
original works. For the other applications of the Wright function, including 
Mikusihski's operational calculus and the theory of integral transforms of 
Hankel type we refer, for example, to Kiryakova [16], Kratzel [17] . Mikusihski 
[27] . and Stankovic [39] . 
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2. Analytical properties 



2.1. Asymptotics 

Probably the most important characteristic of a special function is its asymp- 
totics. In the case of an entire function there are deep relations between its 
asymptotic behaviour in the neighbourhood of its only singular point - the 
essential singularity at z = oo - and other properties of this function, in- 
cluding distribution of its zeros (see, for example, Evgrafov [BJ, Levin |18j). 
It follows from the Stirling asymptotic formula for the gamma function that 
the Wright function 

oo k 

is an entire function of z for p > — 1 and, consequently, as we will see in 
the later parts of our survey, some elements of the general theory of entire 
functions can be applied. 

The complete picture of the asymptotic behaviour of the Wright function 
for large values of z was given by Wright [12] in the case p > and by Wright 
[33j in the case — 1 < p < 0. In both cases he used the method of steepest 
descent and the integral representation 

0(p, p ] z) = 7 L f e C+*C-' C H3 dC) p > _i, p e c (2) 

2m J-Ra 

where Ho denotes the Hankel path in the C-plane with a cut along the 
negative real semi-axis arg£ = n. Formula (J2j) is obtained by substituting 
the Hankel representation for the reciprocal of the gamma function 

! / e c C s d£, seC (3) 



r(s) 2ni Jna 

for s = pk + f3 into (pQ) and changing the order of integration and summation. 

Let us consider at first the case p > 0. 

Theorem 2.1.1 If p > 0, arg(— z) = £, |£| < n, and 

Zl = (plzDV^e^+^+V, Z 2 = {p\z\) l '^e^/ { ' J+1 \ 

then we have 

( / ) (p,p;z) = H(Z 1 ) + H(Z 2 ), (4) 
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where H{Z) is given by 

( M ( 

z m V \z\ M+1 



ff(Z) =^-^ff;(=^= + 0( ' H.^oc (5) 



k m=0 

and the a m , m = 0,1,..., are defined as the coefficients of v 2m in the 
expansion of 

with 

,x f 1 , P + 2 , 09 + 2)00 + 3) 2 H 

In particular, if (3 £ R we get the asymptotic expansion of the Wright 
function c/>(/3, /3; —x) for x — > +oo in the form 

<j>{p,P\ -x) = x p $-®e axPcos * p cos - /?) + ax p sin^) {cj + 0(2T P )} , 

^ (6) 
where p = j^r, & = (1 + p)p 1+p and the constant c\ can be exactly 
evaluated. 

If we exclude from the consideration an arbitrary small angle containing 
the negative real semi-axis, we get a simpler result. 

Theorem 2.1.2 If p > 0, argz = 9, \9\ < n — e, e > 0, and 

Z = {p\z\) 1 /(P+Ve ie /(< ) + 1 \ 

then we have 

<j>(p,ftz)=H(Z), (7) 

where H(z) is given by ([5]). 

In the case p = the Wright function is reduced to the exponential 
function with the constant factor l/T(/3): 

</>(0,P;z)=exp(z)/T(j3), (8) 

which turns out to vanish identically for (3 = —n, n = 0, 1, . . .. 

To formulate the results for the case — 1 < p < we introduce some 
notations. Let 

y = —z, —ir < argz < 7T, —it < argy < tt, (9) 



5 



and let 

K=(l + p)((-p)-'y) 1/(1+p) . (10) 

Theorem 2.1.3 If -1 < p < 0, | argy| < min{§7r(l + p), ir} - e, e > 0, 
then 

4>(p,(3;z) = I(Y), (11) 

where 

I(Y) = Y^e- Y Yl A ^ Y ' m + 0(y- M ) , Y -> oo, (12) 

I m=0 J 

and the coefficients ^4 m , m = 0, 1 . . . are defined by the asymptotic expan- 
sion 

T(l-!3-pt) >jzt (-l) m A m 

27r(-p)-"*(l + p)(i+P)(*+i)r(< + 1) r((l + p)i + /3 + i + m) 

1 \ 



+0 



T{{1 + p)t + l3 + \ + M) 



valid for argt, arg(— pt), and arg(l — p — pt) all lying between — ir and 7r 
and i tending to infinity. 

If —1/3 < p < 0, the only region not covered by Theorem 2.1.3 is the 
neighbourhood of the positive real semi-axis. Here we have the following 
result. 

Theorem 2.1.4 If -1/3 < p < 0, | argz| < ir(l + p) - e, e > 0, then 

<t>(p,(3;z) = I(Y 1 ) + I(Y 2 ), (13) 
where I(Y) is defined by (fT2j) . 

Yi = (1 + p) ((-p)-^e-) 1/(1+P) ,r 2 = (1 + p) ((-p)-^e— , (14) 
hence 

Y\ = Y if — 7r < arg z < 0, and Y2 = Y" if < arg z < n. 

As a consequence we get the asymptotic expansion of the Wright function 
4>(p,P;x) for x — > +00 in the case —1/3 < p < 0, /? G R in the form: 

<Kp,P;x) = x p ^-^e- axPcosnp cos(-K P (- -0)- ax p sin vrp) {c 2 + 0(2T P )} , 

(15) 



6 



where p = j^-, a = (1 + p)(—p) 1+ p and the constant c 2 can be exactly 
evaluated. When —1 < p < —1/3, there is a region of the plane in which 
the expansion is algebraic. 

Theorem 2.1.5 If -1 < p < -1/3, |arg*| < |tt(— 1 - 3p) - e, e > 0, 
then 

Hp, /3; z) = J(z), z^- oo, (16) 

where 

z (p-l-m)/(-p) _ 1 _ M 

J<2) = E (-P)r(m + m + (/3 - m - !)/(-,)) + 0(2 >' 

Finally, the asymptotic expansions of the Wright function in the neigh- 
bourhood of the positive real semi-axis in the case p = —1/3 and in the 
neighbourhood of the lines argz = ±^7r(— 1 — 3p) when — 1 < p < —1/3 are 
given by the following results by Wright. 

Theorem 2.1.6 If p = -1/3, | argz| < tt(1 + p) - e, e > 0, then 

4>(p,l3;z)=I(Y 1 )+I(Y 2 ) + J(z), (18) 

where I(Y) is defined by (JT2D, Y 1 ,Y 2 by CED, and J(z) by (fTTjl . 

Theorem 2.1.7 If -1 < p < -1/3, |argz± ^(-1 - 3p)\ < tt(1 + p) - 
e, e > 0, then 

0(p,/J;*)=I(y) + J(z), (19) 
where I(Y) is defined by CE]) and J(z) by (fT?) . 

The results given above contain the complete description of the asymp- 
totic behaviour of the Wright function for large values of z and for all values 
of the parameters p > — 1, (3 € C. We will use them repeatedly in our 
further discussions. 

2.2. Representations through hypergeometric functions 

Due to the relation 

<Ki,v + i;-\* l ) = (|) "./„(*). (20) 

Wright considered the function 0(/?, /?; z) as a generalization of the Bessel 
function J u (z). In the general case of arbitrary real p > — 1 the Wright 



7 



function is a particular case of the Fox .ff-function (|13|. [2], [16j App. E], 
[381 Chapter 1]): 



H, 



1.0 
0.2 



(0,1),(1-m,p) 



21 



Unfortunately, since the Fox H- function is a very general object this repre- 
sentation is not especially informative. It turns out that if p is a positive 
rational number the Wright function can be represented in terms of the 
more familiar generalized hypergeometric functions. Let be p = n/m with 
positive integers n and m. Substituting s = ms\ into (|21|) and making use 
of the Gauss-Legendre formula for the gamma function 



T(nz) = n n 
we arrive at the representation 



/ N 1 — ^1 

(27T) — 



n— 1 



Y[T(z + -), n = 2.1. 



n 



m 



n — m X 

(27r) 2 m^ri 



2m 



to— 1 

n + 

k=0 



n-1 

n r( 

(=0 



(- z y 



m" l n" 



dsi, 



(22) 



which is equivalent to the representation given by Pathak [30] in terms of 
the Meijer G-function ([161 App. A], [25l Chapter 4]). Here L_oo is a loop 
beginning and ending at — oo, encircling in the positive direction all the 



poles of r(si + — ), k = 0, . . . , m — 1, i.e., the points 
residue theorem and the relation ((2SJ Chapter 3]) 



.A 

m ' 



1 



TO ' 



The 



res. 



<T(z) 



kl 



k = 0,1,2,... 



allow us to represent this integral as a sum of m series of hypergeometric 
type: 



I TO— 1 OO / -i \ 



TO— 1 , 

n r(- 9 -^) 



• rrr(^ + ^ + g + ^- 

11 V n TO ^ n 

i=0 



(-*)' 



m" l n' 



ton 



Using the Gauss-Legendre formula and the recurrence and reflection formu- 
lae for the gamma function 

T(z + l) = zT(z), r(z)r(i-z) 



simrz 
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to simplify the coefficients of the series in the last representation we obtain 
the final formula 

<K-, P\ z) = £ f - Fn +m -i (-; A(n, I + £), A*(m, £±±); , 
m p!l (p + — p) \ n m m m m n n / 

(23) 

where p F 9 ((a) p ; z) is the generalized hypergeometric function ([251 Chap- 
ter 4]), 

1 A; — 1 

A(fc,a) ={ a ,a + a + — —}, A*(k,a) = A(k, a) \ {1}. 

We note that the set A*(k,a) is correctly defined in our case since the 
number 1 is an element of the set A(m, ^—), < p < m — 1. 

The same considerations can be applied in the case of negative rational 
p but under the additional condition that the parameter (3 is also a rational 
number. In particular, we obtain the formulae 

0(-- 5 -n; z) = 1 J Z r(| + n) ^ + n; t --), n = 0, 1, 2, . . . , 

Z 7T Z Z Z 4 

(24) 

1 1 (-l) n 1 1 1 z 2 

0(-2> 2 " n 5 z ) = ^^^2 +n) lFl( 2 +n; 2 ; ~T } ' n = °' lj 2 ' ' ' ' ' (25) 
If n = we get 

^(-io;z) = -^/ 4 , (26) 

ijz) = 4=e- 9/ *. (27) 

Z Z y 7T 

The formula (|26p was given by Stankovic [39]. He also gave the relation 

(x > 0) 

,2 2 1 2 / 4 \ 

rM — , 0; -x s ) = — exp( — ^rjVF i i — — , 

where Wu,v(%) is the Whittaker function satisfying the differential equation 

E^w + (-i + i + £)" r w = a 

The formula (I27h as well as some other particular cases of the Wright func- 
tion with p and f3 rational, — 1 < p < 0, can be found in Mainardi and 
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Tomirotti [23], where a particular case of the Wright function, namely, the 
function 

M(z; (5) = <f>(-P, 1 - P; -z), < f3 < 1 (28) 
has been considered in details. For (3 = 1/q, q = 2, 3, . . . the representation 



1 im -i 



M(z; — ) = - V (-l) n - l Y(n/m) sm( irn/m)F(z; n, m) (29) 

with 



m 7r 

n=l 



F(z;n,m) = ^(-l) fc ( m+1 )(n/m) fc 

k=0 



(mk + n — 1)! 



was given. Here (a)k, k = 0, 1, 2, . . ., denotes the Pochhammer symbol 

(a) fc = r( " + . fc) =a(Q + l)...(Q + fc-l). 
r(a) 

In particular, the formula (|29p gives us for m = 3 the representation 
4>{-\ 2 --,z)=^M(-z/^) 

with the Airy function Ai(z) . Finally, we rewrite the formulae (|24p . ([25 
by using the Kummer formula ( |25[ Chapter 6]) 

1 F 1 (a;c;z) = e z i_Fi(c - a; c; -z) 

in the form 

1 



2 
1 1 



1 ' "l 



n: : , = e ~ z / 4 zP n (z 2 ), n = 0, 1,2, . . . , (30) 

; Z ) = e - z2 /4Q n(z 2 )i n = 0) l,2,..., (31) 



Yy 2' 2 

where P n (z), Qn{z) are polynomials of degree n defined as 

(_l)n+l z 

P n {z) = ^ r(3/2 + n) xF^-n^/2;-), 

IT 4 

Q n {z) = ^r(l/2 + n) iFi(-n; 1/2; J). 

7T 4 
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2.3. Laplace transform pairs related to the Wright function 

In the case p > the Wright function is an entire function of order less than 
1 and consequently its Laplace transform can be obtained by transforming 
term-by-term its Taylor expansion (pTJ) in the origin. As a result we get 
(0 < t < +oo, s£C, < e < \s\, e arbitrarily small) 

/■OO 

<t>(p,P;±t) + C[</>(jf,P;±t);a]= e'^p, 0;±t) dt (32) 

Jo 



*> _ st ~ (±t) k _ ~ (±l) fc r°° _ st k 
o 6 ^ n k\T(pk + f3) dt -^kT(pk + f3)J 6 



k=0 ^ 1 ^> k=0 

where -j- denotes the juxtaposition of a function t/?(i) with its Laplace trans- 
form <f(s), and 

OO fc 

£ -'* ) = £F(^)- a>0 -' 36C - (33) 

is the generalized Mittag-Lemer function. In this case the resulting Laplace 
transform turns out to be analytic, vanishing at infinity and exhibiting an 
essential singularity at s = 0. 

For — 1 < p < the just applied method cannot be used since then 
the Wright function is an entire function of order greater than one. The 
existence of the Laplace transform of the function <p(p, (5; —t), t > 0, follows 
in this case from Theorem 2.1.3, which says us that the function (f>(p,/3;z) 
is exponentially small for large z in a sector of the plane containing the 
negative real semi-axis. To get the transform in this case we use the idea 
given in Mainardi [23]. Recalling the integral representation ([2]) we have 
(-Kp<0) 



OO 



1 



0(p,#-t) + I e- st <t>(p,P;-t)dt= e~ st — e^" C P d( dt 

JO Jo Ml JHa 

= —[ e^-P I" ' e-^ B+ <~ P UtdC (34) 
2iri Jua Jo 

1 f e^C 13 

2m Jko, s + ( p 

again with the generalized Mittag-Leffler function according to (|33p . We use 
here the integral representation (see Djrbashian [2], Gorenflo and Mainardi 
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E ^) = ^~ I T^^C, (35) 



2m Jua ( a ~z 

which is obtained by substituting the Hankel representation (J3j) for the re- 
ciprocal of the gamma function into the series representation fj33|) . 

The relation (|34|) was given in Djrbashian and Bagian [2] (see also Djr- 
bashian [3] ) in the case (5 > as a representation of the generalized Mittag- 
Lefner function in the whole complex plane as a Laplace integral of an entire 
function and without identifying this function as the known Wright function. 
They also gave (in slightly different notations) the more general representa- 
tion 

E a2 ,p 2 {z) = nE aiA (zt a i)t^- 1 <j)(-a2/a 1 ,p 2 -p 1 — ;-t)dt, 
J o «i 

< a 2 < ai, /3i,/? 2 > 0. 

An important particular case of the Laplace transform pair (134p is given by 

M(t;f3)^E f3 (-s), 0</?<l, (36) 

where M(t; (5) is the Mainardi function given by (|28p and 

oo 

E a {z) = E a \(z) = ]T — — , a > (37) 

is the (standard) Mittag-Lefner function. The formula (|36p contains, in 
particular, the well-known Laplace transform pair 

M(t; 1/2) = exp(-t 2 /4) + E 1/2 (-s) = exp(s 2 ) erfc (s), s 6 C. 
Using the relation 

/ t n f(t)dt = ]hn(-l) n — C[f(t);s), 
Jo s^o ds n 

the Laplace transform pair (|34p and the series representation of the gener- 
alized Mittag-Lefner function (|33|) we can compute all the moments of the 
Wright function 4>(p,l3; —t), — 1 < p < in R + : 

/ f>(p, /?; -t) dt = - — -, n G No = {0, 1,2,.. .}. 

For the Mainardi function M(t; (3), < /3 < 1 we obtain from this formula 
the normalization property in R + (n = 0) 
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M(t; /3)dt = l 

i 

and the moments in the form 

poo n \ 

/ t n M(t;P)dt = • neN. 
Jo r(/3n + 1) 

Now we introduce the function (Mainardi |23j ) 

F(z-f3) =4>(-P,0;-z), 0</5< 1, (38) 

which is connected with the function M(z; (3) by the relation 

F(z;P)=PzM(z',P). (39) 
For this function we can prove the relation 

jF(\t~ p \ p) = J^-Af (Ai _/J ; P) -=- exp(-A A < p < 1, A > 0. (40) 



Indeed, following Mainardi [23J and using the integral representation ([2]) we 
get 



£-Mexp(-A^);i] = / e ^ ds = *~» ^ d( (41) 

27TZ 7Ha 27T2 £ JHa 

= lF(At-/ 3 ;/3) = ^M(A^;/3). 

The Laplace transform pair (|40p was formerly given by Pollard [33] and by 
Mikusihski [27]. 

By applying the formula for differentiation of the image of the Laplace 
transform to (|40p we get the Laplace transform pair useful for our further 
discussions: 

pM(At" /3 ; p) -=- s p ~ l exp(-As /3 ), < p < 1, A > 0. (42) 

In the general case, using the same method as in (141 p . we get (see 
Stankovic [39] ) the Laplace transform pair 

^~V0>> # -Ai p ) + s _/3 exp(-As" p ), -1< p < 0, A > 0. 

Stankovic [39] also gave some other Laplace transform pairs related to the 
Wright function including 

i _/3 exp(-t^ p cos(pvr)) sin(/?vr-i _p sin^wr))^^" ft -s p ), -1 < p < 0, /? < 1. 
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2.4. The Wright function as an entire function 
of completely regular growth 

The fact that the function ([1]) is an entire function for all values of the 
parameters p > — 1 and (3 £ C was already known to Wright (Wright [32], 
|44j ) . In the paper Djrbashian and Bagian [JJ (see also Djrbashian [3]) the 
order and type of this function as well as an estimate of its indicator function 
were given for the case — 1 < p < 0. Wright [H] also remarked that the zeros 
of the function (UJ) lie near the positive real semi-axis if— 1/3 < p < and 
near the two lines argz = ±^7r(3/3 + 1) if —1 < p < —1/3. In this paper we 
continue the investigations of the Wright function from the viewpoint of the 
theory of entire functions. We give exact formulae for the order, the type and 
the indicator function of the entire function <fi(p,(3;z) for p > —1, (3 € C. 
On the basis of these results the problem of distribution of zeros of the 
Wright function is considered. In all cases this function is shown to be a 
function of completely regular growth. 

The order and the type of the Wright function are obtained in a straight- 
forward way by using the standard formulae for the order p and the type a 
of an entire function f{z) defined by the power series f(z) = J2k^=o c nZ n 

ralogra , .-,/„ i/_ — 

p = limsup - — - — T — p-, (crep) IF = lim supra " &c n 

n-+oo log(l/|c n |) n^oo 

and the Stirling asymptotic formula 

T(z) = V2^z z ~ l/2 e~ z [l + 0(l/z)}, \ argz\ < vr - e, e > 0, \z\ -» oo. 

We thus obtain the following result. 

Theorem 2.4.1 The Wright function cj>(p, f3; z), p > -1, (3 £ C ((3 ^ 
—n, n = 0, 1, . . . if p = 0) is an entire function of finite order p and the type 
a given by 

p = -J—,a = (l + p)\p\~^. (43) 
1 + p 

Remark 2.4.1 In the case p = the Wright function is reduced to the 
exponential function with the constant factor l/T(/3), which turns out to 

vanish identically for (3 = —n, n = 0,1,.... For all other values of the 

— p 

parameter (3 and p = the formulae (|43p (with a = lim(l + p)\p\ 1+p = 1) 
are still valid. 

The basic characteristic of the growth of an entire function f{z) of finite 
order p in different directions is its indicator function h(9), \8\ < tt defined 
by the equation 
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^)=lim S up 1 ° g|/(r ^ )l . (44) 

To find the indicator function h p (9) of the entire function <p(p, (3; z) of finite 
order p given by (03]) its asymptotics given in Section 2.1 are used. By direct 
evaluations we arrive at the following theorem. 

Theorem 2.4.2 Let p > -1, (3 e C (/? ^ -n, n = 0, 1, . . . if p = 0). 
Then the indicator function h p {6) of the Wright function 4>(p,f3;z) is given 
by the formulae 

h p {6) = acospO, \9\ < vr (45) 

in the case p > 0, 

h p (6) = \- (TCOSP t + e \ ZV/ 0, ( 46 ) 
H { —a cos p(6 — it), < U < it 

in the cases (a) — 1/3 < p < 0, (b) p = —1/2, /3 = — n, n = 0, 1, . . . and (c) 
p = -1/2, /? = 1/2 - n, n = 0, 1, . . ., 

!-cr cos p(7r + 6»), — 7r < 6» < || - 7T, 
0, |f| < 7T- |S, (47) 

—a cos 

in the case — 1 < p < — 1/3 (/3 ^ — n, n = 0, 1, . . . and /3 / 1/2 — n, n = 
0, 1, . . . if p = —1/2), where p and a are the order and type of the Wright 
function, respectively, given by (}4"3"j) . 

Remark 2.4.2 It can be seen from the formulae (f4"5j) . ([4*6]) that the 
indicator function h p (9) of the Wright function cj)(p,(3;z) is reduced to the 
function cos 9 - the indicator function of the exponential function e z - if p — » 
0. This property is not valid for another generalization of the exponential 
function - the Mittag-Leffier function (|37p . Even though 

Ei(z) = e z , 

the indicator function of the Mittag-Leffier function given for < a < 
2, a / 1 by (0 Chapter 2.7]) 



cos8/a, \0\ < ?f, 
0, ^<|#|<^ 



2 

does not coincide with the indicator function of e z if a — > 1. 
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We consider now the problem of distribution of zeros of the Wright func- 
tion in the case p > — 1, (3 G R. To get the asymptotics of zeros of the 
Wright function we use its asymptotic expansions flU), ([6]), (fl"3|) . (fl~5|) . (fT8|) . 
(|19p and the method applied by M.M. Djrbashian in [21 Chapter 1.2] to 
solve the problem of distribution of zeros of the generalized Mittag-Lemer 
function E p u(z). This method consists in finding the asymptotics of zeros 
of the main terms of the asymptotic expansions, applying the Rouche the- 
orem to show that the function under consideration and the main terms of 
its asymptotic expansions have the same number of zeros inside of specially 
chosen contours and after that in estimation of the diameter of the domains 
bounded by the contours. The proofs of the results given below are straight- 
forward but have many technical details and are omitted in this paper. It 
turns out, that in dependence of the value of the parameter p > — 1 and the 
real parameter /3, there are five different situations: 

1) for p > all zeros with large enough absolute values are simple and 
are lying on the negative real semi- axis; 

2) in the case p = the Wright function becomes the exponential func- 
tion with a constant factor (equal to zero if (3 = — n, n = 0, 1, . . .) and it 
has no zeros; 

3) for —1/3 < p < all zeros with large enough absolute values are 
simple and are lying on the positive real semi- axis; 

4) in the cases p = —1/2, (3 = —n, n = 0,1,... and p = —1/2, (3 = 
1/2 — n, n = 0, 1,... the Wright function has exactly 2n + 1 and 2n zeros, 
respectively; 

5) for —1 < p < —1/3 (excluding the case 4)) all zeros with large enough 
absolute values are simple and are lying in the neighbourhoods of the rays 
argz = ±i"7r(— 1 — 3p). 

We now give the precise results. 

Theorem 2.4.3 Let {7fc}i° be the sequence of zeros of the function 
4>(p, (3; z), p > —1/3, p 7^ 0, (3 G R, where \^fk\ < |7fc+i| and each zero 
is counted according to its multiplicity. Then: 

A. In the case p > all zeros with large enough k are simple and are 
lying on the negative real semi-axis. The asymptotic formula 



is true. Here and in the next formulae p and a are the order and type of the 
Wright function given by ()43p , respectively. 




(48) 
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B. In the case — 1/3 < p < all zeros with large enough k are simple, 
lying on the positive real semi-axis and the asymptotic formula 



Ik 



'irk + 7r{p(3- 2=1) 



-a sin irp 



{l + 0(£T 2 )}, A;^+oo (49) 



is true. 

Remark 2.4.3 Combining the representation (|20p with the asymptotic 
formula (08]) we get the known formula (see, for example |40[ p. 506]) for 
asymptotic expansion of the large zeros r/% of the Bessel function J v {z): 

1 1 

r k = ir(k + -v - -) + 0(k^ r ), k — > oo. 

Remark 2.4.4 In the cases p = —1/2, (3 = — n, n = 0, 1, . . . and p = 
— 1/2, (3 = 1/2 — n, n = 0, 1, . . . the Wright function can be represented by 
the formulae ([30]) . (f3Tj) and, consequently, has exactly 2n + 1 and 2n zeros 
in the complex plane, respectively. 

It follows from the asymptotic formulae (jlip . (|16p and (|19p that all 
zeros of the function 0(p, f3; z) in the case —1 < p < —1/3 with large enough 
absolute values are lying inside of the angular domains 



3vr 

iU ' R : ' = ! 71 ~ 2^ 



< e 



where e is any number of the interval (0,min{7r — |^}). Consequently, 
the function <p(p, (3; z) has on the real axis only finitely many zeros. Let 

{ 7 i +) }?° G G« = : 9 > 0}, { 7 jH}f G G<") = {z : 9 < 0} 

be sequences of zeros of the function <f>(p, (3; z) in the upper and lower half- 
plane, respectively, such that |7^ | < 1 7^+i I > |t& I < l7fe+il> anc ^ eacn zero 
is counted according to its multiplicity. 

Theorem 2.4.4 In the case -1 < p < — 1/3 ((3 7^ — n, n = 0, 1, . . . and 
(3 7^ 1/2 — n, n = 0, 1, . . . if p = —1/2) all zeros of the function cp(p, (3; z), (3 £ 
R with large enough k are simple and the asymptotic formula 



% 

is true. 



<*> = e**-® (^) ' {l + O } , fc - +oo (50) 



17 



Summarizing all results concerning the asymptotic behaviour of the 
Wright function, its indicator function and the distribution of its zeros, we 
get the theorem. 

Theorem 2.4.5 The Wright function <fi(p,[3;z), p > — 1 is an entire 
function of completely regular growth. 

We recall ( |18l Chapter 3]) that an entire function f(z) of finite order p 
is called a function of completely regular growth (CRG-function) if for all 
9, \8\ < 7T) there exist a set E$ C R+ and the limit 

lim (51) 

r — i- + oo yP 

res* 



where 



E$ = R + \E e , lim meS ^ n( °' r) = 0. 



It is known ([6j Chapter 2.6]) that zeros of a CRG-function f(z) are regularly 
distributed, namely, they possess the finite angular density 

lim ^ = u(9), (52) 

where n(r, 9) is the number of zeros of f(z) in the sector < arg z < 9, \z\ < 
r and p is the order of f{z). ^,From the other side, the angular density v{9) is 
connected with the indicator function h(9) of a CRG-function. In particular 
(see [U Chapter 2.6]), the jump of h'(6) at 9 = 9q is equal to 2-irpA, where 
A is the density of zeros of f{z) in an arbitrarily small angle containing the 
ray arg z = 9q. 

In our case we get from Theorem 2.4.2, that the derivative of the in- 
dicator function of the Wright function has the jump 2apsimrp at 9 = ir 
for p > 0, the same jump at 9 = for — 1/3 < p < 0, and the jump 
ap at 9 = ±(tt - g) for -1 < p < -1/3 (fi ^ -n, n = 0, 1, . . . and 
(3 7^ 1/2 — n, n = 0, 1, . . . if p = —1/2), where again p and a are the order 
and type of the Wright function, respectively, given by (|4"3l) ; if p = or 
p = —1/2 and either f3 = —n, n = 0, 1, . . ., or (3 = 1/2 — n, n = 0, 1, . . ., the 
derivative of the indicator function has no jumps. As we see, the behaviour 
of the derivative of the indicator function of the Wright function is in accor- 
dance with the distribution of its zeros given by Theorems 2.4.3, 2.4.4 and 
Remark 2.4.4 as predicted by the general theory of the CRG-functions. 
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3. Some applications of the Wright function 



3.1. Asymptotic theory of partitions 

Historically the first application of the Wright function was connected with 
the asymptotic theory of partitions. Extending the results of Hardy and 
Ramanujan about asymptotic expansion of the function p(n), n G N, the 
number of partitions of n, Wright [II] considered the more general problem, 
namely, to find an asymptotic expansion for the function Pk(n), n G N, 
the number of partitions of n into perfect k-th powers. Following Hardy 
and Ramanujan, Wright considered the generating function for the sequence 
{jPfc(l), Pfc(2), • • •} which is given by 



oo 



Then 



fk(z) = Y[(l-z l r l = l + Y,Pk(n)z n , \z\<l. 

1 = 1 71=1 

f \ 1 f fk(z)dz 



the contour C being the periphery of the circle with center in the point z = 
and radius r = 1 — — . Let the contour be divided into a large number of 
small arcs, each associated with a point 

a Ptq = exp(2p7ri/q), p, q G N. 

Taking the arc associated with «o,i = 1 as typical, it can be shown that on 
this arc the generating function fk(z) has the representation 

, , N z j ( l\a /r(l + (l/ib))C(l + (l/Jb))N 

(53) 

where j is a real number depending on k and (,(z) is the Riemann zeta- 
function. Then, on this arc, fk(z) is approximated to by an auxiliary func- 
tion Fk(z), which has a singularity at z = 1 of the type of the right-hand side 
of (|53p . If the z-plane is cut along the interval (1, oo) of the real axis, Ff t (z) 
is regular and one- valued for all values of z except those on the cut. The 
power series for F^{z) has coefficients given in terms of the entire function 
(j)(p, (3; z) and, by using this power series, an asymptotic expansion can be 
found for Pk(n). 

In the paper [41J Wright gave some properties of the function cp(p, ft; z) 
in the case p > 0, including its asymptotics and integral representation ([2]). 
He proved on this base the following two theorems. 
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Theorem 3.1.1 Let a, j3, 76 C, p > 0, m G N, m > ^(7), 

00 

F(z) = a, 0, r , z) := £ (n - 7)^ # a(n - 7 ) p )^- (54) 

n=m 

If a cut is made in the z-plane along the segment (1, 00) of the real axis, then 
F{z) is regular and one-valued in the interior of the region thus defined. 
Theorem 3.1.2 Let 

G(z)=F(z)- X (z), 
where F(z) is defined by (|54p and 

z 7 (a 



(log(l/z))/3 ^V(l0g(lA))V 

If a cut is made in the z-plane along the segment (—00, 0) of the real axis, 
then G{z) is regular and one-valued in the interiour of the region thus de- 
fined. 

We see that the function F(x) has a singularity of the type of x( z ) a t 
2 = 1. In the case of the function F^(z) used to get an asymptotic expansion 
for the function pk(n) the values 

p= i, a = r(i + I)C(i + i),/? = -i 7 -i 

should be taken in the previous two theorems. 
3.2. Fractional diffusion- wave equation 

Another field in which the Wright function plays a very important role is 
that of partial differential equations of fractional order. Following Gorenflo, 
Mainardi and Srivastava [13] and Mainardi |20|-|23| we consider the frac- 
tional diffusion-wave equation which is obtained from the classical diffusion 
or wave equation by replacing the first- or second-order time derivative by 
a fractional derivative of order a with < a < 2: 

d a u(x,t) ^d 2 u(x,t) „ . , 

Here the field variable u = u(x, t) is assumed to be a causal function of time, 
i.e. vanishing for t < 0, and the fractional derivative is taken in the Caputo 
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sense: 

r d n u( x ,t) 

d a u(x,t) = 0t* ' a = nGN, 

{ T(n-a)l {t ~ T) c>r" dr ' 1 < a < n. 

(56) 

We refer to the equation (|55|) as to the fractional diffusion and to the frac- 
tional wave equation in the cases < a < 1 and 1 < a < 2, respectively. The 
difference between these two cases can be seen in the formula for the Laplace 
transform of the Caputo fractional derivative of order a (see Mainardi |23j): 

d a u(x,t) ^ a _ 1 _ k d k u(x,t) 

dfa - s a u(x, s)-}2s Ik \t=o+, n-l<a<n, n G N. 

fe=0 

(57) 

Extending the conventional analysis to the equation (I55h . and denoting 
by a i x ) an d h(x) two given, sufficiently well-behaved functions, the basic 
boundary- value problems can be formulated as follows (0 < a < 1): 

a) Cauchy problem 

u(x, 0+) = 5(2;), -co < x < +oo; n(=Foo,t) = 0, t > 0; (58) 

b) Signalling problem 

u(x,Q+) = 0, x > 0; u(0+,i) = /i(t), u(+oo,t) = 0, i > 0. (59) 

If 1 < a < 2 the initial values of the first time-derivative of the field variable, 
u(x,0+), should be added to to the conditions (|58|) and (f59l) . To ensure the 
continuous dependence of the solutions on the parameter a in the transition 
from a = 1— to a = 1+, we agree to assume u(x, 0+) = 0. 

Since these problems are well studied in the cases a = 1 and a = 2 we 
restrict ourselves in the further considerations to the case 0<a<2, a^l. 
For the sake of convenience we use the abbreviation 

P = f , (60) 

which implies < (3 < 1. 

Let us introduce the Green functions Q c (x,t;/3) and Q s (x, t;/3) for the 
Cauchy and signalling problems for the equation (|55p . respectively, which 
represent the fundamental solutions of these problems (with g(x) = 5{x) in 
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([58]) and h(t) = 6(t) in ([59]) ). Using the Green functions, the solutions of 
the two basic problems can be given, respectively, by 



/+oo 
gc(x-Z>t;0)g(£)dt, (61) 
-oo 

u(x,t;0)= f G s (x,t - r; P)h(r) dr. (62) 
Jo 

To get the Green functions Q c {x,t;f3) and g s (x,t;{3) the technique of the 
Laplace transform is used. We consider at first the Cauchy problem ([58]) for 
the equation ([55]) with g{x) = g c (x,0+; (3) = 5(x) (and G c (x,0+; (3) = 
if 1/2 < (3 < 1). Denoting the Laplace transform of the Green func- 
tion by Q c (x,s; f3) and using the formula (|57|) we arrive after application 
of the Laplace transform to the Cauchy problem { (|55p . (|58p | to the non- 
homogeneous differential equation 

_ s ^Q c = -S(x)s^- 1 , -oo < x < +oo (63) 

with the boundary conditions 

g c (Too,s;(3) = 0. (64) 

The problem { (|63p . ([64p | has a solution (see, for example, Mainardi [23J) 

c {x, s; 0) = 1 — e-CM/^D)**, -oo < * < +oo. (65) 
2v^'s " 

Comparing this relation with the Laplace transform pair (|42p we represent 
the Green function for the Cauchy problem { (|55|) . ([58[) } in the form 

g c (x, t; (3) = —!=-M(r/VV; (3), t > 0, (66) 
2yV\x\ 



where 



\x\t-P 



is the similarity variable and M(z; (3) is the Mainardi function (|28|) given in 
terms of the Wright function. 

For the signalling problem { ([53]) . ([51?]) } (with h(t) = 5(t)) the application 
of the Laplace transform leads to the homogeneous differential equation 

v ^s _ s 2^ s = Q) x > (6?) 

with the boundary conditions 

g s (0+,s-(3) = l, g s (+oo,s;(3) = 0. (68) 
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Solving this equation, we obtain 

G s (x,s;P) = e-Cz/v^)^ x > o. (69) 

Using the Laplace transform pair (|40p we get the Green function Q s (x,t;P) 
for the signalling problem { (|55|) . (|59p | in the form 

g.(x, t; 0) = -^M(r/VV; 0), t > 0, x > 0, (70) 
yut 

where 

r = xt~ fi (71) 
is the similarity variable and M(z;f3) is the Mainardi function (|28p . 

For more results in FPDE we refer, for example, to Engler [5], Fujita [7J, 
Gorenflo and Mainardi p], p2], Mainardi [20]-[23], Podlubny [32], Priiss 
[33], Saichev and Zaslavsky [35], Samko et. al. [36J, Schneider and Wyss 
[37] and by Wyss [35]. 

Some applications of FPDE have been considered in papers by several 
authors including Giona and Roman [9], Hilfer |15j . Mainardi [21], Metzler 
et al. [26], Nigmatullin [28], Pipkin (31], Podlubny [32]. 



3.3. Scale invariant solutions of FPDE 

Let us consider the abstract equation 

F(u) = 0, u = u(x,t). (72) 

First we give some definitions concerning the similarity method. 
Definition 3.3.1 A one-parameter family of scaling transformations, 
denoted by T\, is a transformation of (x, t, u)-space of the form 

x = A a x, i = X b t, u = X c u, (73) 

where a, b, and c are constants and A is a real parameter restricted to an 
open interval / containing A = 1. 

Definition 3.3.2 The equation (|72p is invariant under the one-parameter 
family T\ of scaling transformations (|73p iff T\ takes any solution u of (|72p 
to a solution u of the same equation: 

u = T x u and F(u) = 0. (74) 

Definition 3.3.3 A real-valued function r](x, t, u) is called an invariant 
of the one-parameter family T\, if it is unaffected by the transformations, 
in other words: 

r](T\(x,t,u)) = r](x,t,u) for all X £ I. 
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On the half-space {(x,t,u) : x > 0, t > 0}, the invariants of the family 
of scaling transformations (|T3|) are provided by the functions (see [29 1) 

77! (x, t, u) = xt~ a/b , r] 2 (x, t, u) = r clb u. (75) 

If the equation (|72|) is a second order partial differential equation of the form 

G(x, t, u, u x , u t , u xx , u a , u xt ) = 0, (76) 

and this equation is invariant under T\, given by (fT3|) . then the transforma- 
tion 

u(x,t) = t c/b v(z), z = xt~ a/b (77) 

reduces the equation ([76]) to a second order ordinary differential equation of 
the form 

g(z, v, v', v") = 0. (78) 

For a proof of this fact we refer in the case of general Lie group methods to 
[29j . In some cases it can be easily checked directly. 

Recently, the scale-invariant solutions for the equation (|55|) (with the 
fractional derivative in the Caputo and Riemann-Liouville sense) and for the 
more general time- and space-fractional partial differential equation (with 
the Riemann-Liouville space-fractional derivative of order (3 < 2 instead of 
the second order space derivative in the equation (|55p ) have been obtained 
by Gorenflo, Luchko and Mainardi |13j . Buckwar and Luchko pQ and Luchko 
and Gorenflo [19J, respectively. In all cases these solutions have been given in 
terms of the Wright and the generalized Wright functions. Here we present 
some results from these papers. 

At first we determine a group of scaling transformations for the fractional 
diffusion-wave equation (|55f) on the semi-axis (x > 0) with the Caputo frac- 
tional derivative given by (|56p . We have in this case the following theorem. 

Theorem 3.3.1 Let T\ be a one parameter group of scaling transfor- 
mations for the equation (|55|) of the form T\ o (x,t,u) = (Ax, X b t, X c u). 
Then, 

b = - (79) 
a 

and the invariants of this group T\ are given by the expressions 

m(x,t) = xr l/b = xt~ a/2 , m{x,t,u) = r c,b u = r 7 n (so) 

with a real parameter 7 = ca/2. 
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Remark 3.3.1 We note that the first scale-invariant rji of A80\) coin- 
cides with the similarity variable |7i| ) which was used to define the Green 
function of the signalling boundary-value problem for the equation i t 55\) . It 
is a consequence of the fact that the equation l\55\) is invariant under the 
corresponding group of scaling transformations. 

It follows from the general theory of Lie groups and the previous theorem 
that the scale-invariant solutions of the equation (|55p should have the form 



u(x, t) = t<v{y), y = xt~ a/2 . (81) 

Furthermore, the general theory says that the substitution (|8ip reduces the 
partial integro-differential equation (|55p into an ordinary integro-differential 
equation with the unknown function v{y). 

Theorem 3.3.2 The reduced equation for the scale-invariant solutions 
of the equation (|55p of the form (|8ip is given by 

(*P^ a n+1 ' a v)(y) = Vv"(y), y > 0, (82) 

where the operator in the left-hand side is the Caputo type modification of 
the left-hand sided Erdelyi-Kober fractional differential operator defined for 
< 5, n-l<a<n£Nby 

n— 1 i j 

(*P? a 9)(y) ■■={K^ n - a \{{r + j--u-)g){yl y > 0. (83) 

3=0 

Here 

f ffr fr(« - ly^u-^+^giyu 1 ! 5 ) du, a > 0, 
(K^g)(y) := T ^ h (84) 

[ g{y), a = o 

is the left-hand sided Erdelyi-Kober fractional integral operator. 

Remark 3.3.2 As it follows from the definitions of the Caputo type 
modification of the Erdelyi-Kober fractional differential operator (|83p and 
the Erdelyi-Kober fractional integral operator (|84p in the case a = n G 
N, the equation (|82|) for the scale- invariant solutions is a linear ordinary 
differential equation of order max{n, 2}. In the case a = 1 (the diffusion 
equation) we have 

(*P2' 1 v)(y) = (7 - \y~fa.) v (y) 

and (|82p takes the form 

Vv"{z) + \yv'{y) - jv(y) = 0. (85) 
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In the case a = 2 (the wave equation) we get 



(*P7 1,2 v)(y) = {i - I - y-r){i - y—)v{y) 

ay ay 

= y 2 v"{y) - 2( 7 - l)yv'{y) + 7(7 - l)v(y) 
and (|82p is reduced to the ordinary differential equation of the second order: 

(y 2 - V)v"(y) - 2(7 - l)yv'(y) + 7(7 - l)v{y) = 0. (86) 

The complete discussion of these cases one can find, for example, in [29] . 
The case a = n G N, n > 2 was considered in pp. 

Solving the equation (|83p we get the following theorems. 
Theorem 3.3.3 The scale-invariant solutions of the fractional diffusion 
equation ([55]) (0 < a < 1) have the form 

u{x,t) = C x t^{~l + r -^) (87) 

in the case — 1 < 7, 7 / 0, and 

«(M) = Ci0(~ 1;-JL) + C 2 (88) 

in the case 7 = 0, where y = xt - ^ is the first scale invariant (|80p and Ci, C2 
are arbitrary constants. 

Theorem 3.3.4 The scale-invariant solutions of the fractional wave equa- 
tion (|55|) (1 < a < 2) have the form 



u(x,t) = CW(-|,l + 71--)=) (89) 



■^((-a, 2 -a), (2, 3 + 2- 



in the case 1 — a < 7 < 1, 77^1 — 77^ 0, and 

= Ci0(-|,l;--^) (90) 



/ V~a a v v 2 ~i 2 v 2 \ 

+C 2 f — H~ 2 > 1; ^) " ^-^((-«. 2 - a), (2, 3 - -); IL) J + C 3 
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in the case 7 = 0, where y = xt 2 is the first scale invariant (|80|) . 4>((/j,, a), (v, b); z) 
is the generalized Wright function given by (/z + v > 0) 

<f>{{ji,a),{v,b)\z)\=^ — — ■ — — rr, ft^R, a,kC, (91) 

and Ci, C2, C3 are arbitrary constants. 

For the elements of the theory of the generalized Wright function (19ip 
including its integral representations and asymptotics we refer to Wright 
|43j in the case v > and to Luchko and Gorenflo [19] in the case of one 
of the parameters /J,, v being negative. 

We consider now the equation (|55p on the semi-axis x > with the 
fractional derivative in the Riemann-Liouville sense: 

r d n u(x,t) 

&*u{x,t) _ W 1 > t a = neN, 

^ ~ [ r(n - a) fa ~ T) n ~ a ~ lu (x,r) dr, n - 1< a < n. 

(92) 

Also in this case the scale-invariants of a one parameter group T\ of scaling 
transformations for the equation (|55|) of the form T\o(x, t, u) = (Ax, X b t, \ c u) 
are given by Theorem 3.3.1. 

Following Buckwar and Luchko [1] we restrict ourselves in the further 
discussion to the case of the group T\ of scaling transformations of the 
form T\ o (x,t,u) = (Xx, \ b t, u). Then the scale-invariant solutions of the 
equation (|55p with the Riemann-Liouville fractional derivative (|92p have the 
form 

u(x,t) =v{y), y = xt' a / 2 (93) 



and the substitution (I93p reduces the partial integro-differential equation 
(|55p into an ordinary integro-differential equation with the unknown function 
v(y) given by the following theorem. 

Theorem 3.3.5 The reduced equation for the scale-invariant solutions 
in the form (I93h of the equation (|55p with the Riemann-Liouville fractional 
derivative (|92|1 is given by 

(Py a a ' a v)(y)='Dv"(y), y>0 (94) 

with the left-hand sided Erdelyi-Kober fractional differential operator P$' a 
defined for < 5, n-l<a<n£Nby 
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{P T S ' a 9){y) ■■= (j[ (r +3 ~ Jy^)j (K^ n ~ a g)(y), y > 0. (95) 

Here (Kg' a g)(y) is the left-hand sided Erdelyi-Kober fractional integral op- 
erator (jM|) . 

The solutions of the equation (|94p have been given by Buckwar and 
Luchko pQ for a > 1. 

Theorem 3.3.6 The scale-invariant solutions of the equation (|55p with 
the Riemann-Liouville fractional derivative (|92p in the case 1 < a < 2 have 
the form (y = xt~ a / 2 ) : 

OL i OL i 

u(x,t) = v(y) = C 1 cf ) (--,l,-y/Vv) + C2<P(--,l,y/Vv) (96) 
with arbitrary constants C\, C2. 



Now we consider the case a > 2: 

Theorem 3.3.7 The scale-invariant solutions of the equation (|55p with 
the Riemann-Liouville fractional derivative (|92p in the case a > 2, a N 
have the form (y = xt~ a / 2 ) : 



KM) = «(y) = £<7;y- 2+ - (1+i) 

j=0 



;i,i), (2-1(1+^2) a 

(a - j, a) 



where Cj, < j < [a] are arbitrary constants and p ^ q 
is the generalized Wright function (see [38]): 



(ai,Ai),...,(a p ,Ap) , 
(b 1 ,B 1 )...(b q ,B q ) 
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z 



(ai,Ai),. . . ,(a p ,A p ) 
(b 1 ,B 1 )...(b q ,B q ) ' 



^ nLir(a, + ^) z k 



In the case 2<a = n£Nwe have the following result. 

Theorem 3.3.8 The scale-invariant solutions of the partial differential 
equation (2 < n G N) 

— = t > 0, x > 0, P > 

ot n 
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have the form (y = x/t n ^ 2 ) : 

n-2 

u(x,t) = J2C jy - 2+ ^+i) 2 *i 

3=0 



[1,1), (2-2(1+^,2) a 
(n-j,n) 



+ C n _i 
(99) 



with arbitrary constants CL, < jf < n — 1. 



Finally, following Luchko and Gorenflo [19], we consider the time-and 
space-fractional partial differential equation 

d a u(x,t) &u(x,t) _ . , 

= *>°^>°^> ' ( 10 °) 

where both fractional derivatives are defined in the Riemann-Liouville sense 
& 



Theorem 3.3.9 The invariants of the group T\ of scaling transformations 
under which the equation (|100p is invariant are given by the expressions 

7/1 (x, t, u) = xt~ a/f3 , r? 2 (x, t, u) = t~ J u (101) 
with an arbitrary constant 7. 

Theorem 3.3.10 The transformation 

u(x,i) = fv(y), y = xt~ a/(3 (102) 

reduces the partial differential equation of fractional order (jlOOp to the or- 
dinary differential equation of fractional order of the form 

(PiP'^M = Vy-P{D^v){y), y > 0. (103) 

Here the left-hand sided Erdelyi-Kober fractional differential operator Pj' a 
is given by (|95p and the right-hand sided Erdelyi-Kober fractional differential 
operator D T & '^ is defined for < 5, n — l</3<nGNby 

(Dfg)(y) ■= (f[(r + j + (i T s +p ' n ~ P 9)(y), y > o, (104) 

with the right-hand sided Erdelyi-Kober fractional integral operator 

(r/g)(y):=\ F ^ J ° { } ^ ' (105) 

g{y), = o. 
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Solving the reduced equation we arrive at the following theorem. 

Theorem 3.3.11 Let 



^<a<(3<2, n-1 < P <n£N. 

Then the scale-invariant (according to the transformation (|102p with 7 > 0) 
solutions of the partial differential equation of fractional order (jlOOh have 
the form 

n 

U (x,t) = tiJ2 c i v j(y), v = *t~ a/ ^ (106) 

where 

Vj (y) = y p - j <(>{(-a, 1 + 7 - a + ^j), (J3, 1 + - i); //P), (107) 

the C,-, 1 < j < n are arbitrary real constants, and <fi(((i,a), (v,b);z) is the 
generalized Wright function given by (|9ip . 

Remark 3.3.3 In the case = 2 the scale-invariant solutions of the 
equation (|100p can be expressed in terms of the Wright function. Indeed, let 
us consider the linear combinations of the solutions (jlOTh with y = xt~ a / 2 : 

Ul (x,t) = f'(VVv 1 (y)+v 2 (y))=t^(-^,l+r,y/Vv), 

u 2 (x,t) = t'y{-VVv 1 (y) + v 2 (y)) = t'y ( t>(~,l+ T ,-y/Vv). 
These scale-invariant solutions are given in Theorem 3.3.6 in the case 7 = 0. 

Remark 3.3.4 For < < 1 the equation (jlOOp has only one solu- 
tion which is scale- invariant with respect to the transformation (I102p . This 
solution has the form (y = xt~ a ^) : 

u(x, t) = P Vl (y) = iV _1 0((-a, 1 + 7 - ol + |), (0, 0);y p /V). 



In the case (3=1, this function is expressed in terms of the Wright function 
(y = xt~ a ): 

u(x, t) = fv x {y) = f<f>(-a, 1 + 7; y/V). 
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